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^ Abstract. The new approach to quantum mechanical problems is proposed. Quantum 

; ^ . states are represented in an algebraic program, by lists of variable length, while operators 

^ are well defined functions on these lists. Complete numerical solution of a given system 

^ can then be automatically obtained. The method is applied to Wess-Zumino quantum 

. mechanics and D—2 and D=4 supersymmetric Yang-Mills quantum mechanics with the 

0\ • SU(2) gauge group. Convergence with increasing size of the basis was observed in various 

(N ■ cases. Many old results were confirmed and some new ones, especially for the D=4 system, 

are derived. Preliminary results in higher dimensions are also presented. In particular the 
►..^ spectrum of the zero- volume glueballs in 4 < D < 11 is obtained for the first time. 
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Supersymmetric quantum mechanics attracts a lot of attention for many 
reasons. First, it describes well defined, finite dimensional quantum sys- 
tems which are ideal laboratory to study supersymmetry in various dis- 
guises dH, m. Second, the conjecture of Banks, Fischler, Shenker and 
(D , Susskind triggers additional interest, in particular in the D=10 supersym- 

1^ ' metric Yang-Mills quantum mechanics (SYMQM) as a model of M-theory 

1^. The latter system has been studied by many authors also in lower 
^ \ dimensions and for a variety of gauge groups. It is expected to have a 

;_j ■ continuous spectrum due to the supersymmetry driven cancellations of 

zero energy fluctuations j^. This is in contrast to its non supersymmetric 
version, which for D=4 becomes the zero volume limit of the well known 
pure Yang-Mills gluodynamics 0, §]• Remarkably, in ten (and not less 
than ten) Euclidean dimensions the continuous spectrum of the supersym- 
metric system contains also a localized zero energy state. Existence of such 
a threshold bound state - the supergraviton - is considered as a sine qua 
non condition for the BESS hypothesis [^. Hence establishing it became 
an arena of intensive and ingenious studies |1C, 11]. Interesting numerical 



methods were also developed ||T2| . The large N limit of SYMQM was studied 
in the mean field approximation, and interesting realizations of the black 
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hole thermodynamics were observed [^]. Simplified (i.e. quenched, D=4) 
systems were also studied by lattice Monte Carlo methods for intermediate 



size groups [14|. The asymptotic behavior in was observed and a possible 
evidence for the two phase structure was found. 

Many interesting results have been obtained studying fully reduced (to 
one point) matrix model (for a review see e.g. |l^). In particular, the 
notorious sign problem, which has been hindering Monte Carlo studies may 
have been recently reduced p^ . 

In these lectures I will present a new approach to study quantum me- 
chanical systems p7| ]. The standard hamiltonian formulation of quantum 
mechanics will be implemented in the computer code with the vectors in 
Hilbert space represented by Mathematica lists with a flexible, dynamically 
varying size. Quantum operators become well defined and simple functions 
on these lists. Fermionic degrees of freedom are easily included. In prac- 
tical applications one has to limit the size of the Hilbert space. However 
dependence on the cut-off can be monitored, and the infinite cut-off limit 
extracted in many physically relevant cases. Of course the method becomes 
computationally demanding for larger systems. Until now it has proven 
applicable for up to 27 degrees of freedom. We begin with simple two- 
dimensional examples: Wess-Zumino quantum mechanics and D=2 SYM 
quantum mechanics. Then the new results for yet unsolved D=4 SYMQM 
will be summarized. Finally preliminary findings for higher dimensions, 
including D=10, will be presented. More complete account of most of these 



results is in Ref|17|. 



General hamiltonian methods have been applied before to complete, 
space extended field theories Recently Matsumura and collaborators 
1£] and Pinsky et al. have studied with this technique a variety of partly 



reduced, supersymmetric theories in lower dimensions (see Refs.||2^ and 
references therein). 

2. Quantum mechanics in a PC 

Action of any quantum mechanical observable can be efficiently imple- 
mented in an algebraic program if we use the discrete eigen basis of the 
occupation number operator a'^a 

{\n>}, \n>=^{a^r\0>. (1) 



"n\ 

The bosonic coordinate and momentum operators are 
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and a typical quantum observable can be represented as the multiple actions 
of the basic creation and annihilation operators]^. Fermionic observables 
will be discussed in subsequent Sections. Generalization to more degrees of 
freedom is evident and will be done there. 

Any quantum state is a superposition of arbitrary number, n^, of ele- 
mentary states |n > 

1st >= S7^a/|n(^) >, (3) 
and will be represented as a Mathematica list 

St = {ns, {ai, . . . , a, J, {n^^}, {n(2)}, . . . , {n("^)}}, (4) 

with Us + 2 elements. The first element specifies the number of elementary 
states entering the linear combination, Eq.(^), the second is the sublist 
supplying all complex amplitudes a/, I = l,...,ns, and the remaining 
Hs sublists give the occupation numbers of elementary, basis states. In 
particular, an elementary state \n > is represented by {l,{l},{n}}. 

Next we implement basic operations defined in the Hilbert space: addi- 
tion of two states, multiplication by a number and the scalar product. All 
these can be simply programmed as definite operations on Mathematica 
lists transforming them in accord with the principles of quantum mechanics. 
It is now easy to define the creation and annihilation operators which act 
as a list-valued functions on above lists. This also defines the action of the 
position and momentum operators according to Eq.(|2|). Then we proceed 
to define any quantum observables of interest: hamiltonian, angular mo- 
mentum, generators of gauge transformations, supersymmetry generators, 
etc. 

Now our strategy is clear: given a particular system, define the list 
corresponding to the empty state, then generate a finite basis of Ncut vectors 
and calculate matrix representations of the hamiltonian and other quantum 
operators using above rules. Given that, the complete spectrum and its 
various symmetry properties is obtained by the numerical diagonalization. 

Dependence of our results on the cut-off Ncut can be monitored. In many 
systems studied so far one can extract meaningful (i.e. Ncut = c>o) results 
before the size of the basis becomes unmanageable. 

3. Two two-dimensional systems 

Wess-Zumino quantum mechanics (WZQM) has one complex bosonic vari- 
able (p{t) = x(t)+iy{t) = xi{t) + ix2{t) and two complex Grassmann-valued 
fermions V'«(t), a = 1,2 The hamiltonian reads (with the mass and 



^ The method can be also extended to non polynomial potentials. 
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the coupling set to 1) 

H = \{pi + pD + !</> + + ((1 + 20)^1^2 + h.c.) , (5) 

Bosonic creation and annihilation operators are introduced as in Eq.(^) for 
each (real) degree of freedom. Fermionic ones - fa, fa ~ can be chosen so 
that ipa = fa, ''Pa — fa ™ this casc 0. With these creation operators we 
generate the basis which contains all (here up to two) fermionic quanta and 
maximum Ncut bosonic quanta. The size of such a basis is then 2{Ncut + 
l)(-^c«i+2). With the rules of the previous section one then easily calculates 
matrix representation of the hamiltonian (|5|) in above basis. The spectrum 
obtained by numerical diagonalization shows rather fast restoration of the 
supersymmetry: lowest state tends to zero and higher bosonic and fermionic 
states become degenerate around N^ut^ 10. All this is summarized by the 
Witten index, Fig. 1, which is simply calculated in the energy eigen basis. 
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Figure 1. Witten index for Wess-Zumino quantum mechanics for 4 <Ncut< 15. 

Due to the additional symmetry in the model all states are doubly degen- 
erate, hence the exact value Iw{T) = 2. 

D=2 supersymmetric SU(2) Yang-Mills quantum mechanics. This system 
has gauge invariance and continuous spectrum, both of which pose new 
challenge for our approach. It is described by three real bosonic and three 
fermionic variables Xa{t),'ipa{t),a = 1,2,3. The hamiltonian reads Q 

H = ^PaPa + igeabc'Pixblpc = ^PaPa + gXaGa, (6) 

with the SU(2) gauge generators Ga = (-abc{xhPc — iipl'ipc)- To satisfy the 
gauge invariance we construct the complete basis from gauge invariant 
states. To this end we use four independent gauge invariant creators: (aa) = 




^ For the details of the implementation of fermionic operators see 
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alal., (af) = 4/], (aff) = eahcalfUl, and (///) = eabJlfUt Since the 



'■a "a' 



fermionic number is conserved, c.f. Eq.(|^, the last three operators (and 
the identity) acting on the empty state generate the four "base" states for 
four sectors of the Hilbert space. All higher states are then obtained by the 
recursive action of the (aa) creator only. In this scheme (referred to as the 
4+4+... scheme) the size of the basis is 4:Ncut- Matrix representation of the 
hamiltonian and the spectrum are then automatically obtained as before. 
The spectrum converges towards the supersymmetric multiplets and zero- 
energy ground state, however this time the convergence with the cut-off is 
slower. This is the feature of the continuous spectrum (c.f. next Section). 
Interestingly there exists a scheme of increasing the basis (2+4+4+...) 
where the spectrum of the anticommutator of SUSY generators {Q, Q} 
reveals the exact super symmetry at every finite cut-off [||, ||, 0. bmce 
the latter converges to the hamiltonian at infinite Ncut one is free to declare 
it as the finite cut-off energy operator. 

Witten index vanishes identically due to the paticle-hole symmetry. 
However one can define the Witten index restricted only to the two fermionic 
sectors {F = 0,1 say) which are balanced by super symmetry. Such index 
does not vanish and carries nontrivial information. Very recently M. Cam- 



postrini has calculated the index with much higher cut-off than in [17| 
(see Fig. 2) confirming that it tends to 1/2 at all T ^. If further estab- 
lished, this would provide an example of the continuous spectrum with 
time independent Witten index [10|. 




Figure 2. Restricted Witten index for D=2 SYMQM with cut-offs N^ut < 100. 



I thank M. Campostrini for providing this figure prior to the publication. 
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4. D=4 supersymmetric SU(2) Yang-Mills quantum mechanics 

Dimensionally reduced (in space) system is described by nine bosonic co- 
ordinates xli^{t), i = 1, 2, 3; a = 1, 2, 3 and six independent fermionic coor- 
dinates contained in the Majorana spinor tjj'^{t), a = 1, ...,4. Hamiltonian 
reads |22| 



Hb = ^PaPa + ^^abceadeXlxixlixi + ^-^eabci^lT'^'llJbx'^, (7) 

where is the transpose of the real Majorana spinor, and in D=4 T are 
the Dirac a matrices. We work in the Majorana representation of Ref. [p3| . 

The system has the rotational symmetry, generated by the spin(3) an- 
gular momentum J* = e^^'^ {xlp^ - jtp'^T,^''ipa^ , T,^'' = -j[T^,T''], the 

gauge invariance with the generators Ga = ^ahc {^^bP'c ~ fV'fe'^c) , and N = 

1 supersymmetry with the charges Qa = {T^'4'a)aPa + is^abc{'^'^'^i^a)axlx^. 
Again, expressing bosonic observables in terms of creation and annihilation 
operators is evident. For fermions we use 



/ 



~^'^fa ~ fa ~ fa^ + ^/a^ 

K-ifl-fl + fJ+ifJ J 



(8) 



which gives the quantum hermitean Majorana spinor satisfying {V'a > V'f } = 
^ab ■ Fermionic operators are implemented with the aid of the Jordan- 
Wigner construction using the lexicographic ordering of the double index 
A = {a, a). It turns out that the fermionic number F = fa^ fa is conserved. 
Consequently the hamiltonian can be diagonalized independently in all 
seven fermionic sectors, which is doneQ analogously to the previous cases. 

The spectrum of the theory is shown in FigJ^ A sample of states is 
labeled with their angular momenta j, which are simply assigned using 
the matrix representation of the angular momentum. It is important that 
our cut-off preserves the SO (3) symmetry. In particular the states have 
appropriate degeneracies for each value of j. For current Ncut the spectrum 
extends to Emax ~ 35. 

Sectors with = 6, 5, 4 can be obtained from those with F = 0, 1, 2 by 
the particle-hole transformation. Since our cut-off respects this symmetry 
we indeed observe, exact repetition of the eigenvalues in corresponding 
subspaces. Therefore only first four sectors are displayed in Figj^. Evidently 



* Construction of the basis is similar to the D—2 SYMQM but more involved and will 
not be discussed here. Full details are in Ref.[[l7|. 
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i^jgure 5. Spectrum of D=4 SYMQM. 

the spectrum of D=4 theory is very rich. We begin its discussion with the 
correspondence to the pure Yang-Mihs quantum mechanics. 

The fermionic part of the hamiltonian, Eq.(|^) vanishes on the F = 
sector. Therefore our effective hamiltonian in this sector is nothing but that 
of the zero volume pure Yang- Mills theory considered by Liischer . Indeed 
we observe rather satisfactory convergence of the F = eigen energies to 
the results of Liischer and Miinster M , c.f. Fig.4. 
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Lower energy levels in F = sector compared with the zero volume glueballs 



Another comparison was done together with van Baal who studied F = 
2 sector in a more complex, non compact, supersymmetric model of the 
same family . At the time of writing lO] our results could not be directly 



compared, due to the specific boundary conditions required in [25|. However 
after this conference van Baal has adapted his programs so that they can 
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be applied also to our situation. We have found complete agreement to the 
full Mathematica precision in both F = and F = 2 sectors^ . 
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Figure 5. Cut-off dependence of the lower levels of the D=4 SYMQM in four 
independent fermionic channels {Ncut=^) ■ 



Next we discuss the convergence with the cut-off (c.f. Fig. 5). The first 
striking property is that the Ncut dependence is different in various fermionic 
sectors. For F = 0,1,5, 6 the dependence is weak, hence the lowest states 
have already converged (within some reasonable precision). For F = 2, 3,4 
however, the convergence is slower. We interpret this as an indication that 
the spectrum is discrete in the F = 0, 1, 5, 6 sectors while it is continuous 
in the "fermion/hole-rich" F = 2,3,4 environments. This is based on our 
finding of an interesting and useful correlation between the localizability of 
a state and the rate of the convergence with the cut-off Ncut- It turns out 
that the spectrum of the momentum operator in the cut-off Hilbert space is 
exactly given by the zeroes of the Hermitte polynomials [^] . Consequently 
the convergence of the free spectrum is slow ~ 0{1/Ncut)- On the contrary, 
for the localized (i.e. bound) states the simple argument suggests that the 
rate of convergence is given by the large x asymptotics of the wave function, 
hence it is fast (e.g. exponential). This different nature of the spectra in 
different fermionic sectors provide the reasonable extension of the result of 
Ref.[^ which was mentioned earlier. Since the fermionic modes are crucial 
to generate the continuous spectrum, it is natural that it does not show 
up in the sectors where they do not exist at all, or are largely freezed out 



I thank P. van Baal for providing his results and programs. 
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by Pauli blocking. On the other hand not aU states in F = 2, 3,4 sectors 
have to belong to the continuum. Supersymmetry together with the discrete 
spectrum in = 0, 1 channels implies existence of the normalizable states 
among the continuum of F = 2, 3, 4 states. As one explicit example we 
give the energy of the F = 2,j = 1 state, shown in Fig. 5 (flatter curve 
beginning at E=6). It has definitely weaker dependence on Ncut than the 
others. This indicates that the lowest 12^?,!^ > state is localized. This may 
be a precursor of the more complex phenomenon expected in the D=10 
theory. There, the zero energy localized bound state of Do branes should 
exist at the threshold of the continuous spectrum. Present example suggests 
that one way to distinguish such a state from the continuum may be by the 
different Ncut dependence. Finally our results show that the supersymmetric 
vacuum cannot be in the F = sector. Instead the best candidates are the 
lowest states in the F = 2 and 4 sectors. 
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Figure 6. Supersymmeric images of a sample of eigenstates. 

SUSY on the level of states. Evidently identifying SUSY multiplets 
in Fig.3 looks like a nontrivial task. However present approach allows to 
monitor directly the action of the SUSY generators on the eigenstates of 
the hamiltonian. As any other observables they can be implemented in 
our "quantum algebra" described in the second Section. Alternatively their 
matrix representations are easily obtained once bases are constructed in 
each sector. Then the supersymmetric image of an eigenstate from one 
sector can be decomposed into states from adjacent sectors (the charges 
change fcrmionic number by ±1). This allows identification of the approx- 
imate supersymmetry partners. With the current scheme supersymmetry 
is restored only at infinite Ncut, hence the energies of the SUSY partners 
are not exactly identical. Neither the SUSY images of eigenstates are 100% 
concentrated on single eigenstates. This situation is quantitatively sum- 
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marized in Fig. 6. The lower multiplets can be already identified, while for 
higher states, larger cut-offs are required for better assignment of SUSY 
partners, which is not surprising. 

Witten index. Figure 7 shows the (Euclidean) time dependence of the 
Witten index calculated from our spectrum with up to five bosonic quanta. 
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Figure 7. Witten index of the D=4, supersymmetric Yang-Mills quantum mechanics 
for the number of bosonic quanta bounded by Ncut = 2, 5. The bulk value 1/4 is also 
shown. 



Clearly we are far from the convergence, but some interesting signatures 
are already there. First, one can see how the discontinuity at T = emerges 
even at this early stage. Second, the exponential fall-off common to all 
curves at large (T > 5) times is evident. This is where the lowest state dom- 
inates. Since at this Ncut its energy is not zero, the fall-off is exponential. 
Finally, and most interestingly, we observe the flattening shoulder appear- 
ing at T ~ 2 — 3. This is the signal of the supersymmetric cancellations 
which occur on the average, even though the exact SUSY correspondence 
between individual states does not yet show up. The behavior with Ncut is 
not inconsistent with the exact bulk value 1/4 obtained from the nonabelian 
integrals [10, Clearly higher Ncut are needed, and, what is important, 
they are within the range of present computers. 



5. From D=4 to D=10 space-time dimensions 



To bridge the gap to ten dimensions we study first the simpler, F = case 
in all space dimensions 2 < d = D — 1 < W. Thereby results presented 
here extend for the first time the spectrum of the zero volume glueballs to 
higher dimensions. To this end we use the hamiltonian in Eq.(^) without 
the fermionic part, drop all the fermionic occupation numbers in elementary 
states, Eq.(^), and appropriately extend the ranges of spatial indices. There- 
fore generalization to higher dimension is trivial in our approach, however 
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it puts considerably heavier load on the computer. Hence our results are 
rather preliminary and will be updated soon. Figure 7 shows the energies 
of the first three levels as a function of d. 
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Figure 8. Three lowest glueball states in various dimensions. 



The higher dimensional glueballs are there as expected. The ordering 
of the SO(d) multiplets remains the same as in d=3: scalar, two dimen- 
sional symmetric tensor with dimension g=d(d-|-l)/2-l, and scalar again. 
Our cut-off respects rotational symmetry for all d, hence we observe these 
degeneracies exactly in the spectra. The largest uncertainty (i.e. that of 
the d=9 points) is around 25% at present. Given above range of d one may 
attempt to confront these results with the large d expansions. The eye-ball 
fit suggest the power dependence E ~ dT^'^. 

We conclude that calculations in higher dimensions including d = % 
are perfectly realistic (for SU(2) gauge group) even if much more time 
consuming. Including fermions will increase this demand by approximately 
factor of 16. 



6. Summary 



The new approach to solve quantum mechanical systems on a computer is 
proposed. The Hilbert space of quantum states is implemented algebraically 
in a computer code with states represented as flexible Mathematica lists. 
Any quantum observable is then defined as the list valued function on these 
lists. This allows for the automatic calculation of matrix representations 
and the spectrum. A range of progressively more complicated (also super- 
symmetric) systems from two to ten space-time dimensions was studied 
and many new results were obtained. The method is practically applicable 
to the D=10 supersymmctric SU(2) Yang-Mills quantum mechanics, which 
is our main goal in the future. Extension to few higher gauge groups is 
conceivable but large N require dedicated and new modifications. 



12 



Acknowledgments. I thank Jeff Greensite and Stefan Olejnik for providing 
the opportunity to present this work, for financial support, and for the 
warm atmosphere enjoyed during the Conference. I would like to thank C. 
M. Bender for an instructive discussion which inspired this approach. I also 
thank P. van Baal, P. Breitenlohner and H. Sailer for intensive discussions. 
This work is supported by the Polish Committee for Scientific Research 
under the grants no. PB 2P03B01917 and PB 2P03B09622. 



References 



9. 
10. 

11. 

12. 
13. 



15. 
16. 
17. 
18. 
19. 
20. 
21. 

22. 



23. 
24. 
25. 
26. 



E. Witten, Nucl. Phys. B185/188 (1981) 513. 

F. Cooper, A. Khare and U. Sukhatme, Phys. Rep. 251(1995) 267, |hep-th/940502c . 
M. Claudson and M. B. Halpern, Nucl. Phys. B250 (1985) 689; J. Greensite and 
M. B. Halpern, Nucl. Phys. B242 (1984) 167. 

T. Banks, W. F ishier, S. Shenker and L. Susskind, Phys.Rev. D55 (1997) 6189; 



hep-th/9610043. 



B. de Wit, M. Liischer and H. Nicolai, Nucl. Phys. B320 (1989) 135. 

M. Liischer, Nucl. Phys. B219 (1983) 233. 

M. Liischer and G. Miinster, Nucl. Phys. B232 (1984) 445. 

P. van Baal, Acta Phys. Polon. B20 (1989) 295; also in At the Frontiers of Particle 
Physics - Handbook of QCD, Bor is loffe Festschrif t, vol. 2, ed. M. Shifman (World 
Scientific, Singapore 2001) p.683; [hep-ph/000820e| . 

J. Polchinski, String Theory, Cambridge University Press, Cambridge, 1998. 

S. Sethi and M. Stern, Comm. Math. Phys. 194 (1998) 675; P. Yi, Nucl. Phys. 

B505 (1997) 307; A. V. Smilga, Nucl. Phys. B266 (1986) 45. 

G. Moore, N. Nekrasov and S. Shatashvih, Commun. Math. Phys.209 (2000) 77; 
M.B. Green and M. Gutperle, JHEP 01 (1998) 005. 

W. Krauth and M. Staudacher, Nucl. Phys. B584(2000) 641, |hep-th/0004076 
D. Kabat, G. Lifschytz and D. A. Lowe, Phys. Rev. D64 (2001) 124015 , 



th/0105171 



hep- 



14. R. A. Janik and J. Wosiek, Acta Phys. Polon B32 (2001) 2143, |hep-th/9903121 
P. Bialas an d J. Wosiek, Nucl. Phys. B (Proc. Suppl.) 106 (2002) 968 , |hep- 



lat/01 11034 



H. Aoki et al., Prog. Theor. Phys. Suppl. 134 (1999) 47, hep-th/9908038 



J. Ambjor n, K.N. Anagnos topoulos and A. Krasnitz, JHEP 0106 (2001) 069. 
J. Wosiek, |hep-th/0203"ll6| 



J. Kogut and L. Susskind, Phys. Rev. Dll (1975) 395. 
Y. Matsumura, N. Sakai and T. Sakai, P hys. Rev. D52 (1995) 2446. 



J. R. Hiller, S. Pinsky and U. Trittmann, [hep-th/011215l| ; |hep-th/0106193 , 
M. A. Shifman, ITEP Lectures on Particle Physics and Field Theory, World 
Scientific, Singapore, 1999. 

M. B. Halp ern and C. Schwartz, Int. J. Mod. Phys. A13 (1998) 4367, |he£ 



th/9712133 



C. Itzykson and J.-B. Zuber, Quantum Field Theory, McGraw-Hill, New York, 1980. 
P. Jordan an d E. P. Wigner Z. Phys. 47 (1928) 631. 



P. van Baal, 
J. Trzetrzelewski and J. Wosiek, in preparation. 



[iep-th/0112072. 



